TOEPLITZ AND HANKEL DETERMINANTS WITH SINGULARITIES: 

ANNOUNCEMENT OF RESULTS 



P. DEIFT, A. ITS, AND I. KRASOVSKY 

Abstract. We obtain asymptotics for Toeplitz, Hankel, and Toeplitz+Hankel determi- 
nants whose symbols possess Fisher-Hartwig singularities. Details of the proofs will be 
presented in another publication. 



Let f{z) be a complex- valued function integrable over the unit circle with Fourier coeffi- 
cients 

^27r 

We are interested in the n- dimensional Toeplitz determinant with symbol f{z), 

(1) Dnifi^)) = det(/,_,)^_-io. 

In this paper we present the asymptotics of Dn{f{z)) as n ^ oo in the case when the 
symbol /(e*^) has a fixed number of Fisher-Hartwig singularities [17, 25], i.e., when it has 
the following form on the unit circle: 

m 

(2) f{z) = 6^(^)^^.-0/3. -Q 1^ _ ^^.|2a.^^^_^^(^)^-ft-^ , = ^^e^ Q ^ [0^27r), 



/. = 7^/ f{^'y~^''de, j = 0,±l,±2. 



for some m = 0, 1, . . . , where 



(3) z^=e''\ j = 0,...,m, ^ = % < < ■ ■ ■ < <2t^- 

/ N / N fe'"'^^ < argz < % 

4 g^fi{z)=gp{z) = l 

^ ^ ^ 1 e ""'^J 9j < arg 2; < 27r 

(5) > -1/2, pj eC, j = 0, . . . , m, 

and V{e^^) is a sufficiently smooth function on the unit circle. Here the condition on a 
insures integrability. Note that a single Fisher-Hartwig singularity at zj consists of a root- 
type singularity 

(6) Iz-zA'''^ = 2siJ~^' 



and a jump g/3j{z). A point zj, j = 1, . . . , m is included in (3) if and only if either aj 7^ or 
Pj 7^ (or both); in contrast, the point Zq = 1 is always included (note that gp^^z) = e~*'^^")_ 
Observe that for each j, z^^gp^i^z) is continuous at z = 1, and so for each j each "beta" 
singularity produces a jump only at the point Zj. The factors z- are singled out to simplify 

comparisons with existing literature. Indeed, (2) with the notation h{9) = e^^^'*-* is exactly 
the symbol considered in [17, 2, 3, 4, 5, 6, 7, 10, 11, 12, 15, 16, 30]. We write the symbol, 
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however, in a form with z^i=°^^ factored out. The present way of writing is more natural 
for our analysis. 

On the unit circle V{z) is represented by its Fourier expansion: 

(7) V{z) = J2 Vkz\ Vk = Y Vi^'>~'''d9. 

k=-oo ^ "'O 

The canonical Wiener-Hopf factorization of e^*^^^ is 

First, we address the (essentially known) case when all lie in a single half-closed 
interval of length 1, namely "^Pj G (g — 1/2, g + 1/2], q G M, reproving results of Szego for 
aj = f3j = 0, Widom [30] for (3j = 0, Baser [2] for 3?/?^ = 0, Bottcher and Silbermann [10] 
for \^aj\ < 1/2, \^(3j\ < 1/2, Ehrhardt [16] for \^(3j -^(3k\<l, and other results of these 
authors (see [16] for a review). Note that we write the asymptotics in a form that makes it 
clear which branch of the roots is to be used. 

Theorem 1. Let /(e*^) be defined in (2) and aj ± I3j 7^—1, —2, . . . for j = 0, 1, Then 

as n ^ oo, 



(9) £>„(/) =exp 



j=0 



k=i 

0<j<k<m ^ ^ ^ 



aj(3k-akl3j 



m 



if ^aj > ~, - 3?/3fe| < 1, j. A; = 0, 1, . . . , m, 

where G{x) is Barnes' G-function. 

Remark 2. In the case of a single singularity, i.e., when m = Oorm = l, ao — — 0, the 
theorem implies that the asymptotics (9) hold for 

(10) ^am>~, /3m eC, ± /3„ ^ -1, -2, . . . 

In fact, if there is only one singularity and F = 0, an explicit formula is known [10] for D„(/) 
in terms of the G-f unctions. 

Remark 3. Assume that the function V{z) is analytic. Then the following can be said about 
the remainder term. If all f5j = 0, the error term o(l) = 0{n~^lnn). If there is only one 
singularity the error term is also 0(n~^ Inn). In the general case, the error term depends on 
the differences Pj — Pk- For analytic V{z), our methods would allow us to calculate several 
asymptotic terms rather than just the main one presented in (9) (and also in (22) below). 

Remark 4. If all ^pj G (-1/2,1/2] or all ^pj G [-1/2,1/2), the conditions aj ± pj 7^ 
— 1, — 2, . . . are satisfied automatically as ^aj > —1/2. 
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Remark 5. Since G{—k) = 0, A; = 0, 1, . . . , the formula (9) no longer represents the leading 

asymptotics if aj + f3j or Uj — f3j is a negative integer for some j. A similar situation arises 
in Theorem 9 below if some representations in Jv[ are degenerate. We do not address this 
case in the paper. 

As mentioned above, Theorem 1 was proved by Ehrhardt. We give an independent proof 
of this result using a connection of -Dn(/) with the system of polynomials orthogonal with 
weight f{z) (2) on the unit circle. First, we can show that all Dk{f) 0, k — ko,ko + 1 . . . , 
for some sufficiently large k^. Then the polynomials (f)k{z) — Xk^^ + " " " ? (t>k{z) — Xk^^ + • • • 
of degree k, k — kQ,kQ+ 1, . . . , satisfying 



1 

2^ 



27r 



4>k{z)z ^f{z)de = Xk^^jk, 



1 

2^ 



2n 



)z^f{z)d9 = Xk'Sjk, 



j ^0,l,...,k, 



exist and are given by the following expressions: 



(12) 



(13) 



where 



1 


/oo 

/lO 


foi 

hi 




fok 
flk 




/fc-10 

1 


fk-11 
Z 




/fc-lfc 
^k 


1 


/oo 

/lO 


foi • • • 

/ll ••• 


fok- 
flk- 


-1 1 
-1 z~ 


V^kDk+l 


fkO 


fkl • • • 


fkk- 


-1 z- 



1 f^^ 

fst = ^ I f{z)z-^'-'^de, s,t^ 0,1,..., k. 



We obviously have 



(14) 



Xk 



D 



k+l 



These polynomials satisfy a Riemann-Hilbert problem. We solve the problem asymptotically 
for large n in case of the weight given by (2) with analytic V{z), thus obtaining the large 
n asymptotics of the orthogonal polynomials. The main new feature of the solution is a 
construction of the local parametrix at the points Zj of Fisher-Hartwig singularities. This 
parametrix is given in terms of the confluent hypergeometric function. A study of the 
asymptotic behaviour of the polynomials orthogonal on the unit circle was initiated by 
Szego [28]. Riemann-Hilbert methods developed within the last 20 years allow us to find 
asymptotics of orthogonal polynomials in all regions of the complex plane (see [13] and 
many subsequent works by many authors). Such an analysis of the polynomials with an 
analytic weight on the unit circle was carried out in [26], and for the case of a weight with 
Q!j-singularities but without jumps, in [27]. We provide, therefore, a generalization of these 
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results. Here we present only the following statement we will need below for the analysis of 
determinants. 

Theorem 6. Let /(e*^) be defined in (2), V{z) he analytic in a neighborhood of the unit circle, 
and (f)k{z) = XkZ^ + ' ' ' > (f>k{z) = Xkz'^ + ■ ■ ■ be the corresponding polynomials satisfying (11). 
Assume that |3?/9j — 3fJ/9fc| < 1, aj ± Pj ^ —1, —2, . . . , j, k — 0,1, . . . ,m. Let 



(15) 

Then as n ^ oo, 
(16) xl-i = exp 



5 = maxn2^(^-^^-i). 



2n 



k=0 



i=0 k^j 



^k I 



^3 ^k \ ^k 



Uk r{aj - /3j)r{ak + Pk) b_{zj)b+{zk) 

+0{S^) + 0{6/n)) , 



where 
(17) 



Vj = exp < —ITT 



\P=0 p=j+l 



\7- - 7 



Under the same conditions, 

-2/3,-u„,, r(l + + b+izj) 



(18) </>n(0) = Xn ( J]n 

0=0 



T{aj-(3j) b-{zj) 



+ 



S + 



n 



n 

max 

k n 



(19) 0n(U) Xny2^^n z^ ^^^^^^^^ 6^.,)+^ 

Remark 7. The error terms here are uniform and differentiable in all aj, jSj for jSj in compact 

subsets of the strip 1 3fJ/?j — 3fJ/9fe | < 1, for aj in compact subsets of the half-plane SfJttj > —1/2, 
and outside a neighborhood of the sets ± /9j = —1, —2, .... \i -\- fij — or — /3j — 
for some j, the corresponding terms in the above formulas vanish. 

Remark 8. Note that the terms with n'^(^''~^^~^^ in (16) become larger in absolute value that 
the 1/n term for maxj^k^iPj — Pk) > 1/2. 

Our proof of Theorem 1 uses Theorem 6, similar results for the asymptotics of the or- 
thogonal polynomials and their Cauchy transforms at the points Zj, and a set of differential 
identities for the logarithm of D^, in the spirit of [14, 20, 23]. 

Our next task is to extend the result for arbitrary (5j e C, i.e. for the case when not 
all 3fJ/3j's lie in a single interval of length less than 1. We know from examples (see, e.g. 
[10, 8, 16]) that in general, the formula (1) breaks down. Obviously, the general case can be 
reduced to 3?/3j e (g — l/2,g-|-l/2] by adding integers to Pj. Then, apart from a constant 



n 



n 



max ■ 

k 



n 
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factor, the only change in f{z) is muhiphcation with z^, i E 1^. However, as can be shown, 
the determinants Dn{f{z)) and Dn{z^f{z)) are simply related. For example, for £ = 1, 2, ... , 



(20) 
where 



Dn{z'f{z)) 



1 1.7=1 J- 



$n(0) 
$n(0) 



dz 



$n+l(0) 
Tz^nMO) 



^^n{0) |l^$„+l(0) 



$n+^-l(0) 
^^n+£-l(0) 



l(0) 



_«___<T> 

and ^k{z) = (j)k{z)/xk- Since the Xk, 4>k{0), 0fe(O) for large k are given by Theorem 6, 
and expressions for the derivatives can be found similarly, it is easy to obtain the general 
asymptotic formula for However, this formula is implicit in the sense that one still needs 
to separate the main asymptotic term from the others: e.g., if the dimension i of F„ is larger 
than the number of leading-order terms in (18), the obvious candidate for the leading order 
in Fn vanishes (this is not the case in the simplest situation given by Theorem 14). We 
outline below how we resolve this problem. 

Following [8, 16], define a so-called representation of a symbol. Namely, for f{z) given by 
(2) replace f3j by f3j + nj, rij G Z if Zj is a singularity (i.e., if either f3j ^ or aj ^ or 
both: we set no = if 2:0 = 1 is not a singularity). The integers rij are arbitrary subject to 
the condition Y^jLo'^^j ~ 0- ^ slightly different notation from [8, 16], we call the resulting 
function f{z; no, . . . , Um) a representation of f{z). (The original f{z) is also a representation 
corresponding to no = • • • = rim — 0.) Obviously, all representations of f{z) differ only by 
multiplicative constants. We have 



(21) 



f{z) = Y[ z]' X f{z; no, ... , n„). 



3=0 



We are interested in the representations (characterized by {nj)Y=o) of / such that ^^o(^'^i + 
rijY is minimal. There is a finite number of such representations and we provide an algorithm 

for finding them explicitly (see Remark 11). We call the set of them Ai. Furthermore, we 
call a representation degenerate if aj + {j3j + rij) or aj — {j3j + rij) is a negative integer for 
some j. We call Ai non- degenerate if it contains no degenerate representations. We prove 



Theorem 9. Let f{z) be given in (2), ^aj > -1/2, j3j G C, j 
non- degenerate. Then, as n ^ 00, 



0, 1, . . . , m. Let M be 



(22) 



^»(/)=E n 



TZ{f{z; no, ... , n„,))(l + o(l)). 



where the sum is over all representations in M.. Each TZ{f{z;no, . . . ,nm)) stands for the 
right-hand side of the formula (9), without the error term, corresponding to f{z; no, . . . , n^). 

Remark 10. This theorem was conjectured by Basor and Tracy [8]. The case when the 
representation minimizing YlJ=o{^Pj + ^j)^ unique, i.e. there is only one term in the 
sum (22), was proved by Ehrhardt [16]. Note that this case happens if and only if there 
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exist such rij that ^f3j + Uj belong to a half-open interval of length 1 for all j = 0, . . . , m: 
see next Remark. Thus, Theorem 9 in this case follows from Theorem 1 applied to this 
representation. 

Remark 11. The set A4 can be characterized as follows. Suppose that the seminorm = 
maxj,fc|SR/3j - ^(3k\ > 1. Then, writing = + 1, (3^^^ = A - 1, and (3^^^ = (3j if 
j ^ s,t, where (3s is one of the beta-parameters with 3?/5s = minj^Pj, (3t is one of the 
beta-parameters with 3?/3t = maXj'SiPj, we see that ||/?*-^-'|| < and / corresponding to 
P^^^ is a representation. After a finite number, say r, of such transformations we reduce an 
arbitrary set of Pj to the situation for which either < 1 or = 1. Note that 

further transformations do not change the seminorm in the second case, while in the first 
case the seminorm oscillates periodically taking 2 values, WP^''^ and 2 — Thus all 

the symbols of type (2) belong to 2 distinct classes: the first, for which < 1, and the 

second, for which \\(3^^^\\ — 1. For symbols of the first class, M has only one member with 
beta-parameters (3''^K Indeed, writing bj = ^[3j, if —1/2 < b^J^^ — g < 1/2 for some g e R and 
all j, then for any (%)^o ^^ch that Yl]Lo % — and not all kj are zero, we have 
(23) 

m mm m m m m 

j=0 j=0 j=0 ,7=0 j=0 j=0 j=0 

where the first inequality is strict as at least one kj > 0. For symbols of the second class, we 
can find g G M such that — 1/2 < b^p —q < 1/2 for all j. Equation (23) in this case holds with 
">" sign replaced by ">". Clearly, there are several representations in M. in this case (they 
correspond to the equalities in (23)) and adding 1 to one of I3^s^ with = miuj b'^p = q-1/2 

while subtracting 1 from one of (3[^^ with 6^^''^ = max^ b^ = g + 1/2 provides the way to find 
all of them. 

A simple explicit sufficient, but obviously not necessary, condition for M. to have only one 
member is that all 3fJ/3j mod 1 be different. 

Remark 12. The situation when all Uj ± (3j are nonnegative integers, which was considered 
by Bottcher and Silbermann in [11], is a particular case of the above theorem. 

Remark 13. The case when all the representations of / are degenerate (not only those in 
M.) was considered by Ehrhardt [16] who found that in this case Dn{f) = 0(e"^'n'^'), where 
r is any real number. We can reproduce this result by our methods but do not present it 
here. 

(r) 

We prove Theorem 9 in the following way. Consider the set (3j constructed in Remark 
11. We have to consider only the second class, i.e. = 1. We then have, relabelhng 

(r) 

(3j according to increasing real part, 

= ■ ■ = < < ■ ■< < = • • • = 

for some p,i > 0. Here m' = m + 1 if 2; = 1 is a singularity, otherwise m' = m. Now consider 
the symbol (not a representation of /) / of type (2) with beta-parameters denoted by (3 and 
given by Pj = Pj"^ for j = 1, . . . ,m' — i, and Pj = Pj'^^ — 1 for j = m' — £+1, . . . , m'. It is easy 
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to see that the original symbol / has [^^^) representations in M. obtained by shifting any i 

out of i+p parameters with the smallest real part to the right by 1. Thus, f{z) — cz^ f{z), 
where c is a simple constant factor. To find the asymptotics of Dn{f), we now use (20) with 
/ replaced by /. The l.h.s. is then c~^Dn{f). In the r.h.s. we have a factor Dn{f) to 
which Theorem 1 is applicable since \\f3\\ < 1, and an i x i determinant F„ involving the 
polynomials orthogonal w.r.t. / (for simplicity, consider V{z) analytic in a neighborhood of 
the unit circle) . It is a crucial fact that the size £ of this determinant is less than the number 
of terms, i + p, in the expansion of (t>n{^)IXn of the same largest order 0{n~'^^^^~^) (see (18) 
with (5j replaced by This fact enables us to extract the leading asymptotic contribution 
to Fn (resolving the problem mentioned above). The asymptotics of F„ and D^if) combine 
together and produce (22). 

We will now discuss a simple particular case of Theorem 9 and present a direct independent 
proof in this case. 

Theorem 14 (A particular case of Theorem 9). Let the symbol /^(^) be obtained from 
f{z) (2) by replacing one with f3j^ ± 1 for some fixed < jo < m. Let 3?Q;j > 
e (-1/2, 1/2], i = 0, 1, . . . ,m. Then 

(24) Dr.{f^{z)) = z7"^I?„(/(z)), D^{r{z)) = zl^D^{f{z)). 

All Xn 

These formulas together with (18,19,16,9) yield the following asymptotic description of Dn{f^) 
Let there be more than one singular points Zj and all aj + pj 0. For /^(-z), let jSj^, 
p — 1, . . . , s be such that they have the same real part which is strictly less than the real 
parts of all the other [3j, i.e. 3?/3j^ — ■ ■ ■ = ^Pj^ < minjyj^^...j^ ^Pj. For f~{z) let one Pj^, 
p = 1, . . . ,s be such that = ■ ■ ■ — ^Pj^ > maxj^j^^,„j^ ^Pj. Then the asymptotics of 
Dn{f^) are given by the following: 

s s 

(25) D^(f^) = zj; Yl + D^(n = 4 E + 

p=i p=i 

where 'R'j,± is the right-hand side of (9) (without the error term) in which Pj is replaced by 
Pj ±1, respectively. 

Proof. For simplicity, we present the proof only for V{z) analytic in a neighborhood of the 
unit circle. Consider the case of f~{z). It corresponds to one of the Pj shifted inside the 
interval (-3/2,-1/2]. Since 

we see that 

ri.z) = -Zj^z-^f{z). 

Therefore, using the identity (an analogue of (20) for £ — —1) 

(26) D^{z-'f{z)) = i-ir^DM, 

Xn 
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we obtain 

Dn{f-{Z)) = {-Z,XDn{z-'f{z)) = zl^D^U{z)). 

An 

If, for some ji, j2, • • • , js, we have that = • • • = 'SiPj^ > maxj^j^^ . j^ ^Pj, then we see 
from (19) that only the addends with n^^-'i"^, . . . , n^^^""^ give contributions to the main 
asymptotic term of Dn{f~{z)). Using the relation G{1 + x) = r{x)G{x), we obtain the 
formula (25) for Dn{f~{z)). The case of f'^{z) is similar. □ 

Example 15. In [8] Basor and Tracy noticed a simple example of a symbol of type (2) for 
which the asymptotics of the determinant can be computed directly, but are very different 
from (9). Up to a constant, the symbol is 



(27) f{e^') 



-i, <9 <Tr 
i, TT < 6* < 27r 



We can represent / as a symbol with /3-singularities Pq = 1/2, (3i = —1/2 at the points 
Zq — 1 and Zi — —1, respectively: 

(28) fiz) = g,,^/2{z)g-^,-^/2iz)e^-^' 

We see that f{z) = f~{z) and jo — 1. Therefore by the first part of Theorem 14, we have 

Dnifiz)) = (-l)"M^D„(/(^)), 
Xn 

where (t)n{z), Xn, Dn{f{z)) correspond to f{z) given by (2) with m — 1, zq — 1, zi — e^'^, 
pQ = Pi = 1/2, ao = «! = 0. 

Substituting (19,9) into the above equation for Dn{f{z)), we obtain 

l + (-l)" /2^., .^^2^,o/oA2, 



(29) DMiz)) = \ ' y-G(l/2)^G(3/2)^(l + o(l)), 

which is the answer found in [8]. 

Alternatively, noting that s = 2, ji = jo = 1 and j2 = and using (25) we obtain 

DMiz)) ^ {-ln{-lrn^,^+no,-). 

Since 7^l _ = Uo- = {2n)-^/'^G{l/2fG{3/2f{l + o(l)), we obtain the same result. 

As noted by Basor and Tracy, f{z) has a different representation of type (2), namely, with 
Po = —1/2, Pi = 1/2, and we can write 

(30) f(z) = -^i,_v2(^)^-i,i/2(^)e-'^/^. 

This fact was the origin of their conjecture. In the notation of Theorem 9, the symbol (28) 
has the two representations minimizing ^^j^oii^J + %)^) with no — rii — and the other 
with no — —1, ni — 1. 

Note that in the case Yl%=o(^j = we can always assume that ^Pj G [—1/2,1/2]. The 
beta-singularities then are just piece-wise constant (step-like) functions. This case is relevant 
for our next result, which is on Hankel determinants. 



TOEPLITZ DETERMINANTS 9 

Let w{x) be an integrable complex- valued function on the interval [—1, 1]. Then the Hankel 
determinant with symbol w{x) is 

n-l 



j,k=0 



(31) Dn{w{x)) = det (^j x^+^w{x)dx^ 
Define w{x) for a fixed r = 0, 1, . . . as follows: 

r+l 

(32) w(x) = e^(") n 1^ - "^a;,-(x) 

j=0 

e— 5Rx>a'' 3^/5. e (-1/2,1/2], 



/3o = /5r+i = 0, 3?aj- > --, j = 0,1,..., r + l. 



1 

'2' 

where U{x) is a sufficiently smooth function on the interval [—1, 1]. Note that wc set f]^ = 
(3r+i = without loss of generality as the functions ujq, uJr+i are just constants on (—1, 1). 
We prove 

Theorem 16. Let w{x) be defined in (32). Then as n ^ oo, 
(33) Dn{w) = D„(i)e[("+'^°+"'-+i)^o-°°^(i)-°-+i^(-i)+5 Er=i kv,^] 

r 

X 4-(^"+°o+"r+i+Eo<j<fc<r+i aj«fc+E;=i /3|) (•27r)"o+"'-+in^(°o+°r+i)+Ej"=i(a --/^p 



X n 1^^- - A,r^("^-'=+^^-^'=) |a,.A, - 1 + ^(1-A2)(l-Ai) 

0<j<A;<r+l 

''g(1 + 2«o)G(1 + 2«,+O^P ' G(l + 2a,) + ^W) ' 

^+1 ]^ / 1 1\ 

A = ^afe, 3?a,>--, 3fJ/3,- e , j = 0, 1, . . . , r+l, /J^ = = 0, 

fe=o ^ ^ 

where F(e*^) = (cos ^^), Zj = e*^^, Aj = cos 6'^, j = 0, . . . , r + 1, and the functions b±{z) are 
defined in (8). 

Remark 17. Dn{l) is an explicitly computable determinant related to the Legendre polyno- 
mials (it can also be written as a Selberg integral), c.f. [31], 

Remark 18. Since Pj enter the symbol only via e^^'"^^ , the theorem describes the general case 
with the exception of the situation when some 3?/3j = 1/2 mod 1. 
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To prove Theorem 16 we use the fact that w{x) can be generated by a particular class of 
functions f{z) given by (2). Namely, we can find an even function f of 9 (/(e*^) = 
6 e [0,27r)) such that 

(35) w(x) = f^-V, x^cose, xe\-l,l]. 

Ismfc'l 

It turns out that we must have m — 2r + 1, 9o — 0, 9r+i — tt, 9m+i-j — 27T — 9j, j — 1, . . .r. If 
we denote the beta-parameters of f{z) by Pj, we obtain (3q — /S^+i = 0, /3j = — Z^^+i-j = —Pj, 
j = 1, . . . ,r. In particular, Yl^=o = as remarked above. 

We obtain Theorem 16 from Theorem 1 and asymptotics for the orthogonal polynomials on 
the unit circle with weight f{z) using the following connection we establish between Hankel 
and Toeplitz determinants: 

(36) D4H.)r ^ "=i-2.-. 

where w{x) and f{z) are related by (35). 

Remark 19. Asymptotics for a subset of symbols (32) which satisfy a symmetry condition 
and have a certain behaviour at the end-points ±1 were found by Basor and Ehrhardt in [4]. 
They use relations between Hankel and Toeplitz determinants which are less general than 

(36) but do not involve polynomials. For some other related results, see [19, 24]. 

Remark 20. Asymptotics of a Hankel determinant when some (or all) of f3j have the real 
part 1/2 can be easily obtained. For the corresponding f{z) this implies that certain = 
-1/2 and ^Pm+i-j = 1/2 and the rest 3?3fe e (-1/2, 1/2). Thus, Theorem 9 can be used 
to estimate D2n{f{z)). For the asymptotics of (f)2n{z) in this case we need an additional 
"correction" term which is now 0{n~^^^~^) — 0{1). 

Remark 21. One can obtain the asymptotics of the polynomials orthogonal on the interval 
[—1, 1] with weight (32) by using our results for the polynomials (pki^) orthogonal with the 
corresponding even weight on the unit circle and a Szego relation which maps the latter 
polynomials to the former ones. 

Our final task is to present asymptotics for the so-called Toeplitz+Hankel determinants. 
We consider the four most important ones appearing in the theory of classical groups and its 
applications to random matrices and statistical mechanics (see, e.g., [1, 18, 22]) defined in 
terms of the Fourier coefficients of an even / (evenness implies the matrices are symmetric) 
as follows: 

(37) det{fj_k + fj+k)]^Lo, det(/j_fc - /,+fc+2)2fcio' ^^^ifj-k ^ fj+k+iYj^Lo- 
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There are simple relations [29, 21, 1] between the determinants (37) and Hankel determinants 
on [—1, 1] with added singularities at the end-points, namely, 

„2 



2n''-2n+2 



(38) det(/,_. + f.+kYrklo = —^DMe''^''^)/VT^), 

(39) det(/,_, - fj+k+2)rkLo = — i^„(/(e^'(^))A/r^), 



,2 



(40) det(/,_, + f,^,-,i)r,Lo = /(e^'(^))' '^^'^ 



1 — X 



„2 



(41) det(/,_. - -^io = (^/(^^'^^Vl^ 

where / is even, and x = cos 9. It is easily seen that if f{z) is an (even) function of type 
(2) then the corresponding symbols of Hankel determinants belong to the class (32). Thus 
a combination of the above formulas and Theorem 16 gives the following 

Theorem 22. Let f{z) be defined in (2) with the condition /(e*^) = /(e~*^). Let 9r+i = tt. 
Then as n ^ oo, 

(42) D'^'^^ = e"^o + |[(Qo+ar-+i+s+t)Vb-(ao+^)V(l)-K+i+t)V(-l)+Er=i'=^fe''] 

r 

^ 2(l-«-*)"+P+E^=iK-/5|)-^(°o+ar+i+s+t)^+5(ao+ar+i+s+t)^^(ag+a^ 

l<j<A;<r 

r _ 



JJ 1 1 - I \1-Zj\ -2"^^ ("0+^) 1 1 + 1 -2"^- ("^+1+*) 



X 



^i(ao+a.+i+.+t+i)G(i/2)2 ^ G(l + a,+/3,)G(l + a,-/3,) ^ . . 

11 r7n^9«.^ K^^o{L)), 



G{1 + ao + s)G{l + ar+i + t) jj^ G{1 + 2aj) 

~ 1 

A = -(cKo + Qir+l + S + t) + ^ Qfj, 





3?aj 


1 

>-2' 




0' 


= 0,l,...,r + l. 




Here 














(43) 




det(/,-_ 






-2n + 2, s = i = 


1 

~2 


(44) 


Dl'-'' = 


det(/,_ 


-k - fj+k+2yik=0^ 


with p 


= 0, . = t = l 




(45) 




det(/,-_ 


-k ± /j+fe+l)j,feio' 


with p 


1 


-4 
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Remark 23. For the case 3?/3j = 1/2 see Remark 20 above. 

Remark 24. For the determinant det(/j_fc + fj+k+i)]k=o when the symbol has no 

a singularities at z = ±1 and \^Pj\ < 1/2, the asymptotics were obtained in [5] (see also [6] 
if / is non-even, aj = 0). Note that for symbols without singularities, i.e. for f(z) — e^^^\ 
the asymptotics of all the above Toeplitz+Hankel determinants (and related more general 
ones) were found recently in [7]. 
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